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1. INTRODUCTION 


The effects of atmospheric turbulence in both horizontal and 
near-horizontal flight, during the return of the Space Shuttle, are important 
for determining design, control, and "pi lot-in-the-loop" effects. A non- 
recursive model (based on von Karman spectra) for atmospheric turbulence along 
the flight path of the Shuttle Orbiter has been developed which provides for 
simulation of instantaneous vertical and horizontal gusts at the vehicle center- 
of-gravity, and also for simulation of instantaneous gust gradients. Based on 
this model the time series for both gusts and gust gradients have been generated 
and stored on a series of magnetic tapes which are entitled Snuttle Simulation 
Turbulence Tapes (SSTT). The time series are designed to represent atmospheric 
turbulence from ground level to an altitude of 10,000 meters. 

The purpose of this document is to provide any potential user of 
the SSTT with an appropriate description of the characteristics of the simu- 
lated turbulence stored on the tapes, as well as instructions regarding their 
proper use. Section 2 contains a discussion of the characteristics of the 
turbulence series, including the spectral shape, cutoff frequencies, and 
variation of turbulence parameters with altitude. Information regarding 
the tapes and their use is presented in Section 3. References cited are 
Included in Section 4, Appendices A and B present the results of spectral 
and statistical analyses of the SSTT. 
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2. CHARACTERISTICS OF SIMULATED TURBULENCE 


The non-recursivc turbulence node! used to generate the SSTT is 
based on von Kantian spectra with finite upper limits corresponding to the 
dimensions of the Space Shuttle, relative to the scale of turbulence in the 
atmosphere. Because the scale of turoulence increases with altitude while 
the dimensions of the Space Shuttle are fixed, the finite upper limits of 
the vor. Karman spectra increase with altitude. In order L o take into account 
these spectral changes, for each gust or gust gr.>. iient there are actually 
four time series corresponding to four altitude bands extending from ground 
level to 10,0'JO meters, as indicated in Table 2-1. A r.ore detailed description 
of the characteristics of the turbulence is provided in the subsections which 
follow. 


TABLE 2-1. TURBULENCE PARAMETERS FOR 
ALTITUDE BANDS 


Band 

# 

lower 

Limit 

(m) 

Upper 

Limit 

(m) 

Tire Interval 
(dimensionless) 

Finite Limit 
of Spectrum 
(dimensionless) 

Von Karman Length 
Scalp (m) 

— 

T 

‘1MAX 

f! 2MAX 

r ‘3MAX 


L 2 

r 

L 3 

1 

0 

30 

.6018 

5.22 

3.38 

7.22 



18 

r 

30 

100 

.2300 

13.66 

11.14 

31.27 


99 

78 


100 

762 

.09431 

j 

33.31 

33.76 

120.27 

300 



n 

762 

10000 

.05309 

59.18 

59.97 

213.68 

533 

533 

Li2iJ 


2.1 TURBULENCE GENERATION PROCEDURE 

The six types of SSTT are presented in Table 2-2. For each gust and 
gust gradient series indicated in the table, the generation procedure involved 
convolving a discrete white noise signal of unit variance with a discrete ap- 
proximation of the Impulse response function corresponding to the appropriate, 
dimensionless spectrum (11 . Each of the resulting series consists of 8500 
discrete signals. The time interval associated with each series was based on 
the maximum frequency for which the simulation procedure is cons dered valid. 
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TAtNLt 2-2. lYits C'r MfllL AltO lUKbULtNlt 


Type 

Curresponding 

Spectrum 

Conn entu 

U 1 

: 11 

longitudinal gust 

u 2 

*2? 

transverse gust 

U 3 

'33 

vertical r .t 

jUg/ 

A 

'22/11 

yaw 

)u^/jx^ 

‘'31/11 

pi ten 

:<u 3 /4x 2 

1 31-22 

rol 1 


these time intervals and the corresponding limiting frequencies, are 
included in Table 2-1. The limiting frequencies were calculated according 
to the relation 


imax 


a V*i 


( 2 - 1 ) 


when* 

a - l.U‘> 

- scale of the i-th component of atmosphori< turbulence 
^ = characteristic length of Space Shuttle in i-th direction 

Values of for the four bands are given in Table 2-1 while the characteristic 

lengths, are presented in Table 2-3. It is important to note that the 

values of correspond to von Karman spectra as opposed to Dryden spectra and 

were derived from the local isotropy requirement discussed In subsection 
3 . 2 . 



TABIF 2-3. CHARACTERISTIC DIMFNS IONS 
OF THE SPACE SHUFFLE 


Characteri Stic 
Length 

Mjjni tude 
..(ft]_ (il) 

Explanation 

h 

39 . 56* 12.06* 

mean aerodynamic chord 

e 2 

39.05 11.9 

1/2 wingspan 

*3 

10.95 3.34 

1/2 fuse’age thickness 

2.2 DIMENSIONLESS El 

J-RGr CONTENT 



The total dinensionless enerqy content of each tine series for 
each altitude band was established by integrating the corresponding spectra 
over the appropriate finite limits indicatr-d in Table 2*1. The resulting energy 
content is presented in Table 2-4. As might be expected the total dimensionless 
energy content of each of the turbulent gust series is less than unity. The 

TABLE 2-4. DIMENSIONLESS ENERGY CONTENT 
FOR GUSTS AND GUST GRADIENTS 


Spectrum 

Altitude Band 

1 

2 

3 

4 

♦ii 

.5388 

.7841 

.8956 

.9298 

*22 

.5772 

.7942 

.8952 

.9296 

*33 

.5225 

.7646 

.8809 

.9197 

*22/11 

1.2832 

6.6484 

24.768 

54.125 

*33/11 

1.1321 

5.9699 

22.644 

49.528 

*33/22 
i 

.7049 

L_ . 

4.9954 

22.893 

50.057 


This value Is slightly greater than that previously reported (1] and is 
derived from a more accurate source £2) • 
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dimension 1 sss energ. * content for each gust gtadtent, however, is not limited 
in such a manner and range as high as 54. 125. For both gusts and gust gradients 
the total energy content increase with altitude because of similar increases 
in the limits of integration. 

2.3 VALIDATION OF SIMULATED TURBULENCE 

A soectral analysis of ea .h of t H e dimensionless time series has been 
carried out by means of a Fast foiiier Transfer’ F; T4 13! . The resilts, which are 
presented in Appendix A, demonstrate that the simulated turbule -e possesses 
the proper von Kansan spectral .haracter i sties . 

All of the dimensionless tire series nave *l»o been analy:c'l static 
tically to determine the gust and gust gradient proDability density functions. 

As shown in Appendix r . t w e results of tr.o.e ar ilyse^ indicate *'rt o**h fne 
simulated gusts and gust gradients are normally distributed, with ne ir-zens 
means and standard deviations consistent with the energy content presented 
in Table 2-4. 


Actually the term 


"energy" 


Is not precise when dealing with gust gradients. 
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3. USE OF SIMULATED TURBUl f.NCE TAPES 


The dimensionless simulated turbulence time series are stored on 
six magnetic tapes as surrinar » zed in Table 3-1. Each tape actually contains 
four time series corresponding to the four altitude bands described in Section 
2. The appropriate procedures for reading the tapes are presented in sub- 
section 3.1, while the proper method for converting the time series from 
dimensionless to dimensional form is described in subsection 3.2. 


TABLE 3-1. INDEX OF SHUTTLE SIMULATED 
TURBULENCE TAPES (SSTT) 


Tape 

Time Spries 

Comments 

SSTT - 1 

Uj - gust 

longitudinal gust 

SSTT-2 

U£ - gust 

transverse gust 

SSTT-3 

u^ - gust 

vertical gust 

SSTT-4 

- gust gradient 

yaw 

SSTT-5 

3u^/3Xj - gust gradient 

pi tch 

SSTT-6 

3uj/3x 2 - gust gradient 

roll 


3.1 READING THE TIME SERIES TAPES 

The four time series on each tape are stored In parallel In 
4-word logical records and are correlated (i.e., at any point in the time 
series the 4 turbulence values are all generated from the same string of 
random numbers). Each time series consists of 8500 elements and thus each 
tape contains 8500 4-word records. Pertinent characteristics of the tapes 
are summarized In Table 3-2. 

The first record on each tape contains a 36-character alphanumeric 
descriptor. The second record contains the time series identification number 
(1-6), the number of points in the time series, and the dimensionless genera- 
tion time step size for each altitude band. The format for this record Is 
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T A 1 


\U\ 


ric TAPt 


JI c T r t 

I M W * v A i Ml< 


N.rtn*r o f t»\icls - 
Ho.V*Oi tv;:*' 
char a ' * * t t v i't . 

R • J.ng .if sity 


9 

‘vV • i It-’- i 

7-bit a-*:; 

r.j . sits 1..1 i ti h 


"2I1J,4(IX,M4.7)". Fellow:'-; those two rt-corS, rhe ti«- sen. , is u ' 
in 4-word re ords a-, or*. . ir ,siy d*'Hi iKJ and in trie ton;. at "4*l»,LIl. 

The order of storage in ee h record is from lowest to niy 'St altit.i_> 

Thus the first word ;n each record corresponds to bar - ! -I, the Stvond to b.r .• 

*2, etc. 

In the J 'Im! use of lie t i ■ o series tays the sarplina f r-_ . y 

may be different fro;: the freguer , at which the tapes wc geroraf *J. i* 

fact, the dv .ns ionics s sampling fre;*ency ..ill generally be i.r lo. Ino> 

fore it will be necessary to interpolate the t v e series in ord-.- to v_t valu-s 
at the proper points in dimensional time. Either zero-order or Mrst-order 
interpolation should be used. Also, as time progresses and altitude changes, 
it will be necessary to switch altitude bands in the tiro series consistent 
with Table 2-1. Because of the manner in which the 4 time series are generated, 
any discontinuity due to switching tire series should be minimal. 


3.2 CONVERSION TO DIMENSIONAL VAUjES 

The dimensionless time series on each tape must be converted to 
dimensional form before actual use in a simulation exercise. The conversion 
process generally involves multiplication and/or division by thp appropriate 
turbulence parameters. For dimensionless gusts, u^, the corresponding standard 
deviation, Oj, should be used. Thus 

Uj * OjUj (3-1 

where 

* 

Uj * dimensional gust 
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For dimensionless gust gradient, — — , the parameters o. and L. are used. 

ax. i j 

Thus J 


3u 


3 * j 


!i 

V x j 


( 3 - 2 ) 


where 




dimensional gust gradient 


Conversion of dimensionless time step, T, involves vehicle velocity, V, and the 
x-compone.it of turbulence scale, L^. Thus 

At* = aL T/V ( 3 ‘ 3 ) 

where 

★ 

At = dimensional time steo 


It is important to note that because both L. and V vary with altitude, 

* A 

the dimensional time step At is not a constant. To obtain dimensional 

* 

ti.ie, t , a summation process is involved as follows: 


where 


f * 

L it n 
n=0 n 

N 

aTL L 1n /V n 
n=0 ln n 


■In 


MV 

v(z n ) 

altitude at nth step 


( 3 - 4 ) 



The variation of the turbulence standard deviation, ' , with altituae is 
presented in Table 3-3. The sa-o table contains the turbulence scale, 1 ^ , 
as a function of altitude. Notice should be taken that the values for and 
presented are designed for use with von Karman spectral models and therefore 
differ somewhat from previously tabulated values ( 4 ] which we re designed for 
use with Dryden spectra. The von Karman o. and L. have been computed based on 
the requirement for local isotropy, which ran be expressed as 

2 2 2 

_°1 _ 4 

I 2/3 " : 2/3 " 7 2/3 (3 ' 5 ) 

L 1 4 L 3 

This method of computation is consistent with the established procedure [ 5 ]. 

The vehicle speed, V, is a function of altitude but also nay vary 
from one trajectory to another. Table 3-4 provides values of 

V as a function of altitude. 



TABLE 3-3. VARIATION OF VON KARMAN STANDARD DEVIATION 
AND LENGTH SCALE WITH ALTITUDE 








TABlf 3-4. VARIATION OF SHUT TLt VELOCITY 
WITH AlTITUni (6| 


— ~ ‘ ■ 

* 

ALTITUDE 

V 

(»*) 

(m/^ec) 

100 

^ lx > 

TOO 

1 66 

SOU 

1W 

;ooo 

170 

■ym 

18H 

cu'i 10 

200 

BOO*' 

240 

lOOOO 

300 


4. 
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APPENDIX A 


SPECTRAL ANALYSIS OF SIMULATED TURBULENCE 

By means of a Fast Fourier Transform [3] spectral analyses of all 
simulated turbulence have been performed. The results are presented in 
dimensionless form in Figures A-l through A-24. Table A-l provides a summary 
of these figures. Also included in each figure is the theoretical von Karman 
spectra. The agreement between the theoretical spectra and the computed spectra 
is quite satisfactory. 


TABLE A-l. 

MATRIX OF 

SPECTRAL 

ANALYSIS 

FIGURES 

SERIES 

TYPE 

ALTITUDE BAND 

1 

2 

3 

4 

U 1 

A-l 

A-2 

A- 3 

A- 4 

U 2 

A- 5 

A-6 

A- 7 

A-8 

u 3 

A-9 

A- 10 

A- 11 

A- 12 

3u2/SXj 

A- 13 

A- 14 

A- 15 

A-16 

3u^/3x j 

A- 17 

A-18 

A- 19 

A-20 

auyaxg 

A-21 

A-22 

A-23 

A-24 


A-l 




Gust Spectrum, Altitude Band #1 Figure A-2. u, - Gust Spectrum, Altitude Band • 2 






Figure A-3. u, - Gust Spectrum, Altitude Rand #3 Figure A-4. u. - Gust Spectrum, Altitude Band #4 
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Figure A-5. u-, - Gust Spectrum, Altitude 8and #1 Figure A-6. u 2 - Gust Spectrum, Altitude Band #2 
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Figure A-7. u- - Gust Spectrum, Altitude Band <H Figure A-8. u £ - Gust Spectrum, Altitude Band #4 
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Figure A-13. - Gust Gradient Spectrum , Figure A-14. - Gust Gradient Spectrum 

Altitude Band *1 Altitude Band #2 
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Figure A-19. 3uj/*x. - Gust Grac’ent Spectra-., figure A- 20. ':u Grader? jpect r, j 

Altitude Band ¥3 Air tide *A 
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Figure A-21. 3u 3 /5x 2 - Gust Gradient Spectrum. Figure A-22. 3u 3 /Sx 2 - Gust Gradient Spectrum 

Altitude Band *1 Altitude Band * 2 
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APPENDIX B 


STATISTICAL ANALYSIS OF SIMULATED TURBULENCE 

By means of standard statistical analysis procedures each of the 
SSTT has been analyzed to determine its mean value, standard deviation, and 
probability density distribution. The resulting mean values are pres<-rted 
in Table 9-1 while Table B-2 contains the resulting standard deviations. 

As expected all mean values were near zero. The standard deviations 
represent the square root of the energy content. The ratio of the square 
rojt ji each theoreticd energy content (from Table 2-4) to the corresponding 
standard deviation (from Table B-2) is presented in Table B-3. The agreement 
appears quite satisfactory. 

The gust and gust gradient probability density distributions are 
presented in Figures B-l through B-24 in accordance with Table 8-4. In each 
figure the corresponding normal distribution is also presented. The results 
indicate that both the gust and gust gradient time series are very close to 
normal distributions. 


TABLE B-l. MEAN VALUE OF GUST AND GUST GRADIENTS 


SERIES 

ALTITUDE BAND 

TYPE 

1 

2 

3 

4 

U 1 

-.006109 

-.009637 

-.015505 

-.020793 

U 2 

-.005858 

-.010306 

-.015392, 

-.018100 

u 3 

-.005597 

-.010238 

-.015364 • 

-.018100 


-.000015 

-.000269 

-.002197 

-.006183 

3Uj/3Xj 

-.000014 

-.00027 

-.002198 

-.006188 

3u 3 /3x 2 

-006585 

- 7019388 

-.043538 

. ... { £ v , i -S 

-.064540 




TABLE B-2. STANDARD DEVIATION OF GUST AND GUST GRADIENTS 


SERIES 

TYPE 

ALTITUDE BA*.n 

1 

2 

3 

4 

U i 

. 724352 

.868642 

.927039 

.9466/2 

u 2 

.756057 

.834176 

.936133 

.94/373 

u 3 

.719277 

.867329 

.928437 

.942133 

Ttu^/'Xj 

1.13989/ 

2.591114 

4. 99 J053 

7.387333 

v v i 

1.670499 

2.45426 

4.777666 

7.063241 

' u 3 / '*2 

„ ^ - J 

.836263 

2.225^3., 

4.766939 

7.049476 


TABLE B-3. RATIO OF SQUARE R<" : C, TH'ORtUCAL 

ENERGY CONTENT* TO THE OBSERVED STANDARD 
DEVIATION 1- 


SERIES 

TYPE 

ALTITUDE BAND 

1 

2 

3 

4 

U l 

1.0134 

1.0194 

1.0208 

1.0186 

u 2 

1.0048 

1.0079 

1.0107 

1.0177 

u 3 

l.OOSO 

1.0082 

1.0109 

1.0179 


.9938 

.9951 

.9955 

.9959 

*V a *i 

.9939 

.9955 

.9960 

.9964 

auj/axj 

1.0039 

1.0041 

1.0037 

1.0036 


y~ ~ — — - - ■ 

Theoretical energy content token from Table 2-4. 

^Observed standard deviation taken from Table B-2. 
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B-5 


Figure B-2. u, - Gust Probability Density Distribution, Altitude Band 




Figure B-3. u, - Gust Probability Density Distribution, Altitude Banc »3 
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Figure B-4. u, « Gust Probability Density Distribution, Altitude Band 





Figure B-6. 




Figure B-7. u ? - Gust Probability Density ’J' r t»*i button , Altitude Band 



Figure 8-8. u 9 - Gust Probability Density Distribution, Altitude Band M 







Figure B-10. u.. - Gust Probability Density Di stri bution , Altitude Band ti? 






Figure B-ll. - Gust Probability Density Distribution, Altitude Band *3 



Figure B-l 2. - Gust Probability berisity D’ stribufor- , Alti t jCo 







Figure B-14. 3u~/5x. - Gust Gradient Pro'-.'Lr ' ’ ty '< ity 3 • sf’Djticr 
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B-18 


Figure B- 15 . bujbx. - Gust Gradient Probability Density Distribution, Altitude 
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Figure 8-16. ju ? /:.x. - Gust Gradient Pror.abi 1 . :v .'er sit/ ;tr : yjtion 




Figure B-17. au,/3x, - Gust Gradient p ro!,di>nUy Density 01-;tr1i*u*. on, A’t'tude Band 
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Figure B- 18. ju-j/3x, - Gust Gradient I'robciL' 1 1 :y Dent’ 
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Figure 8-19. - Cost Gra^i*.»nt IV 1 i ry Dtr^'ty 3' str ii^t 'or , Altitude 
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Figure 3-21. »u »/*x 0 - 'Just , ii , ai_" i trNt ; r» » * / L-i stf v^t ’ . Altitude 
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Figure B-22. au^/3x ? - Gust Gradient Probabi ’ i :y density Z' stributior , A'titude Sand 
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Figure B-23. - Gust Gradient Pr ubaL-' 1 < ‘.y Density 1 • s* r : ; t' r , A’titude Banc 
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